Non-Boltzmann behaviour in models of interacting neutrinos 
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We reconsider the question of the relative importance of single particle effects and correlations 
in the solvable interacting neutrino models introduced by Friedland and Lunardini and by Bell, 
Rawlinson and Sawyer. We show, by an exact calculation, that the two particle correlations are not 
"small", and that they dominate the time evolution in these models, in spite of indications to the 
contrary from the rate of equilibration. This result holds even after the model in generalized from 
the original 2 flavor case to A'^ flavors. The failure of the Boltzmann single particle approximation in 
this model is tentatively attributed to the simplicity of the model, in particular to the assumption 
that all neutrinos in the initial state are in flavor eigenstates. 

PACS numbers: 14.60Pq, 05.30.-d 



I. INTRODUCTION 

This paper is an extension of the work of BMcK and 
10 in collaboration with Alexander Friedland j|. Prelim- 
inary accounts of some of the results have appeared in 
the PhD thesis of 10 [2] and the BSc honours report of 
JQi. Details of the calculations omitted here, can be 
found in the thesis and the report. An analysis of the 
two flavor case was briefly reported earlier 0, 

In the early 1990s a number of authors, including 
Thomson and McKellar[i, 0, S, i, E, HII, Pantelone[i2 
and Sigl and Rafelt [l4| . derived neutrino kinetic 
equations which took into account neutrino mixing and 
interactions. They were a set of coupled equations involv- 
ing the number densities (or the single particle density 
matrices) for the different species of neutrinos and anti- 
neutrinos. These kinetic equations were derived from the 
quantum statistical mechanical equations for the many 
particle density matrix by a process similar to the deriva- 
tion of the quantum Boltzmannn equation. In particluar 
it was assumed that many body correlations do not play a 
significant role in the time evolution of the single particle 
density matrix, although, as we discuss below, this is not 
the only assumption in the derivation of the neutrino ki- 
netic equations (or the quantum Boltzmannn equation). 
Neutrino kinetic equations have proved to be very useful 
in studies of the influence of neutrinos on the evolution 
of the early universe [H, [l6l [TtI . [isj . and are now in use 
in the study of supernovae[l9j|. 

In 2003, Friedland and Lunardini [20|. and Bell, Rawl- 
inson and Sawyer [23| re-examined the validity of the 
single-particle approximation inherent in these kinetic 
equations. Part of the motivation was a concern that 
the entanglement induced by the interactions could make 
correlations so important that the single particle approx- 
imation could be invalidated. They studied a simpli- 



fied model which allowed exact solutions of the many 
body problem, and examined the rate of equilibration 
of the system. This discussion was extended by Fried- 
land and Lunardini [22| and by Friedland, McKellar and 
Okuniewiczfl|. The crucial point in the discussion is the 
observation that the incoherent equilibration time and 
the coherent equilibration times scale with the number 
of particles in very different ways: 
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Here 5 is a suitably normalized interaction strength. 
In a large system coherent equilibration is much faster 
than incoherent equilibration, and thus the equilibration 
timescale of the many body system can be used as a 
proxy to determine the validity or otherwise of the sin- 
gle particle approximation, in as much as one does not 
have coherent equilibration in such an approximation. 
Bell, Rawlinson and Sawyer found a coherent timescale in 
their numerical modeling, but Friedland and Lunardini, 
and Friedland, McKellar and Okuniewicz found incoher- 
ent timescales. While these three papers studied systems 
that differed in detail, they lead to the conclusion that, 
at least in many circumstances, the single particle ap- 
proximation, and thus the neutrino kinetic equations or 
the quantum Boltzmann equation, are valid. 

In this paper we demonstrate that this conclusion is 
false. The fiaw in the logic is that while coherent equi- 
libration timescales are an indication that correlations 
are important, incoherent equilbration timescales do not 
necessarily imply that correlations are insignificant. 

Because the model of Friedland, McKellar and Oku- 
niewicz is exactly solvable, we can explicitly calculate the 
two body correlations in this model, and we find that, in 
the case when the number of fiavors, A^, is two, they 
are significant. Moreover we verify that the time evolu- 
tion of the system is determined completely and correctly 
by the two body correlations through the BBGKY equa- 
tions. We go on to demonstrate that the one particle 
approximation, i.e. the quantum Boltzmannn equation, 
with the initial condition that all particles are in flavor 
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eigenstates, leads to the erroneous conclusion that the 
single particle density matrix is constant in time. We 
show that this final conclusion remains valid in the M 
flavor generalization of the model. 

In section II we review the model of Friedland, McKel- 
lar and Okunicwicz, including the M flavor generaliza- 
tion to SU (Af) and we then organize the paper along the 
above lines — section III shows that the two body corre- 
lations are significant, section IV demonstrates that the 
two body correlations and the BBGKY equations cor- 
rectly determine the time evolution of the one particle 
density matrix, and section V shows how the Boltzmannn 
equation fails to determine the time evolution of the sin- 
gle particle density matrix. Finally in section VI we re- 
view our results and discuss their implications. 



II. THE MODEL OF FRIEDLAND, MCKELLAR 
AND OKUNIEWICZ 

Friedland, McKellar and Okuniewicz used the model of 
Friedland and Lunardini, and generalized it to obtain the 
solution to a system of neutrinos with unequal numbers 
of flavors. So that the system could be easily solved ex- 
plicitly and exactly it was restricted to two neutrino types 
and the spatial or momentum degrees of freedom were ig- 
nored. We extend the model to more than two neutrino 
flavors, but the SU{Af) analogue of the 6-j coefficients of 
SU(2) are not so well studied, making the explicit analy- 
sis more difficult to complete. Nevertheless, some general 
results can be obtained for the Af flavor generalization. 
The advantage of the model is that is allows an exact 
investigation of the dynamics of flavor conversion. 

With two neutrino flavors we could regard the neutrino 
system as equivalent to the up-down quark system (or the 
proton neutron system) and use the language of isospin in 
talking about the system. While the the earlier papers 
used the language of spin, because we are considering 
the many flavor generalisation from time to time in this 
paper, we will use the flavor symmetry language, and for 
Af = 2, the language of isospin. 

Neglecting the spatial degrees of freedom (and thus 
the Dirac matrices), the weak interaction of a pair of 
neutrinos inducing flavor mixing in the Af flavor case is 

Hw ^ g {SlkSxp. + S,fj,SxK) , (3) 

where g is a coupling constant. This Hamiltonian acts 
on a product of two A/'-dimensional vectors, the neutrino 
wavefunctions in a flavor-space representation. We in- 
troduce the A/" X TV Gell-Mann matrices [2l],[2i] XjJ e 
{0, 1, . . . , A/"^ — 1}, with the normalization 

Tr(A,Afe) = 2S,k (4) 

and the conventions that 



Ao = VWIa^, TrAj = for j ^ 0, (5) 



where 1;^^ is the Af x Af unit matrix. The diagonal ma- 
trices are, for n e {2, . . . ,A/} 



A„2_i = J— — diag{l,l,...,l,-(n- 1),0, ...,0} 

(6) 

with n — 1 ones and Af — n zeros along the diagonal. For 
the non-diagonal matrices we choose a notation which 
simply extends the Gell-Mann notation for Af — 3. For 
n e {1, . . . ,Af}, and r e {1, . . . , n}, 

{^n^-l+2r-l}i:j — ^i.rSj.n+l + Si^n+lSj,n (7) 
{Ki^-l+2r}i,j — {^i,rSj,n+l — ^i,n+lSj,r) ■ (8) 

Using the identity 

^im'^Ak = \ji:S,Jx^, + - ^ (Aj),,«(Aj)a^ 1 (9) 

the pair Hamiltonian can be rewritten as 

Hw(..) = .,(^ + ii^) (10) 
= 4(F..F,f- '^-y + " ].(n) 

The "F-spin", F = A/2, is the SU{Af) equivalent of 
isospin. The interaction Hamiltonian for N neutrinos is 
then 

N-l N . . s 
1=1 j=i+l ^ ^ 

which is in turn a linear function of the square of the 
total F-spin, 

N 

F = Ef.' (13) 

i=i 

the possible values of which are completely determined by 
the initial configuration of the neutrino flavors. As is 
a constant of the motion, the interaction Hamiltonian is 
completely determined by the eigenstates and eigenvalues 
of the system, which are well defined. 

In the particular case Af = 2, the eigenvalues are com- 
pletely determined by the the total isospin T and the 
number of particles N, 

ET = [gT{T+l)+g^N{N-2)]. (14) 

The eigenstates '^T,M,a are not uniquely determined by 
T, M. One also needs to specify M, the projection of 
the total isospin on the quantization axis, and additional 
quantum numbers a, which specify the different ways 
that states of total isospin T can be constructed. Note 
that, because we are ignoring the position (or momen- 
tum) variables, we are unable to generate states which 
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obey Fermi statistics — any two isospin up particles 
are necessarily identical and thus antisymmetrizing them 
produces a null state, as required by the Pauli principle. 
Take as an initial state a state with definite numbers of 
isospin up and isospin down particles, labelled by U and 
D respectively, 

\-^{t = 0)) = \TuMu) ® \TdMd) = \TuTd;MuMd) ■ 

(15) 

Because we know the eigenstates and the eigenvalues 
of the Hamiltonian of Eq. , we can easily write down 
how this state evolves with time by expanding the initial 
state in terms of basis states with definite total isospin 
and its projection, (T, M), which thus have a well defined 
time dependence: 

\^{t)) = ^e-'*s(^) (r,M| MuMd;TuTd) 

T (16) 
\TuTd]TM). 

To get the one body density matrix for particle 1, which 
we denote as p\^^ [1^, we need to single out that par- 
ticular particle. To do this we decouple it from the rest, 
writing for the initial 

\TuTd) ^\{hku)Tu,TD) (17) 

where Ii — 1/2 is the isospin of the one particle, ku ~ 
T(7 — 1/2 is the total angular momentum of the rest of 
the spin-up particles. 

Note that we have specified that the particle we are 
singling out is initially spin-up. Of course, an analogous 
derivation can be done for a particle that is initially spin- 
down. 

To construct the density matrix of the total system, we 
use the basis mim) in which ku and Td are coupled 
to isospin /, with projection m, but to get the time de- 
pendence we must relate this to states with a fixed total 
isospin T. The result is 

Pit) = |*(t)) {^{t)\ 

- ^ cxp{-ttg[T{T + 1) - T'(T' + 1)]} 

T,T' Jj\lij[ 

(2/„ + 1)^(2/ +l)(2/' + l) 

{TuTd;TM\ MuMd) {TuTd;T'M\ MuMd) 

r 1/2 Tu - 1/2 Tu\{ 1/2 Tu - 1/2 Tu \ 

\td t I ]\Td r r j 

mim) m'im'\ . 

(18) 

To obtain the one-body density matrix, we take the trace 



over the other particles to get 

p['\t)^ exp{-^^5[T(^+l)-T'(T' + l)]} 

T.T\lJij[ 
m.mi ,m'^ 

(2/i + 1)(2/+ 1) 

{TuTd;TM\ MuMd) {TuTd;T'M\ MuMd) 
{IiI;TM\mim) {hI\T'M\ m'^m) 

r 1/2 Tu - 1/2 Tu\{ 1/2 Tu - 1/2 Tu \ 
\Td T I }\Td T' I I 

|/imi) (/im'il . 

(19) 

Eq. pgp is for a particle that initially has spin up. For a 
particle that initially has spin down, one need only inter- 
change the subscripts U with D. This one body density 
matrix can now be used to calculate the probability that 
the particle has spin up at time or the expectation 
value of the spin at time t. It was the basic tool used 
by Friedland, McKellar and Okuniewicz to calculate the 
time evolution of single particle probabilities. 

We wish to calculate two body correlation functions, 
and also test the quantum BBGKY and Boltzmannn 
equations. For these purposes we also require the two 
body density matrix, which is just as straightforward to 
calculate. 

If the two particles initially have the same isospin pro- 
jection, the derivation follows the above, with the differ- 
ence being that one must decouple two spin-up particles 
instead of one. The final result is that one need only 
make the following changes to Eq. p^ . 

• replace Ii = 1/2 with I12 = 1 where I12 is the 
summation of the angular momentum of the two 
particles i.e. I12 — h + h (because both particles 
have isospin up, the pair must have isospin 1), 

• in the 6j symbols replace Tj/ — 1/2 with Tj/ — 1 , 

• replace mi with mi2 where TO12 = -^Ii2, —I12 + 
1, I12 — 1, /12 is the isospin projection of the pair 
of particles, 

in order to obtain the two body density matrix. When 
both isospins are initially up the result is 

P^uHt)- E exp{-ztg[r(T+l)-T'(r' + l)]} 

m,mi2 ,'>TT'i2 

(2/12 + 1) (2/ +1) 

{TuTd;TM\MuMd) {TuTd]T' M\MuMd) 

(/i2/;TM|mi2m) (/12/; T'M\m\2m) 

r 1 Tu-l Tu\{ 1 Tu-l Tu\ 

\td t I ]\Td r I j 

\Ii2m12) (/l2"l'l2l • 

(20) 
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For particles that initially have opposite spins the anal- 
ysis is a little different. In this case, decouple one particle 
from the isospin up system, and one from the isospin- 
down system. The Njj — 1 remaining isospin up particles 
have a total isospin kjj — Tjj — 1/2, which couples with 
ji to give the isospin Tu. The isospin down particles are 
similarly recoupled. Then the single particle isospins Ii 
and I2 are coupled to an isospin I12 , kjj and kjj are cou- 
pled to the isospin /, and I12 and / are coupled to T. 
Standard angular momentum manipulations then allow 
us to construct the two body density matrix for this case, 
which now involves 9-j symbols: 

Pi2{t)= exp{-ztg[r(T+l)-T'(T' + l)]} 

T,T',/,7i2,/(2 

m,mi2 ,m.j^2 

(2Tt/ + l)(2Tz3 + l)(2/+l) 

X [(2/12 + 1) (2/(2 + 1)] 5 

{TuTd;TM\MuMd) (TuTd; T'M\MuMd) 
(/12/; mi2m\TM) m\^m\T' M) 

h ku Tu ) ( h ku Tu 1 
h ko Td > < h ko To > 

I12 I T j [ /12 / r J 

\I12m12) {I'i2m[2\ ■ 

(21) 

With these basic ingredients we can proceed to obtain 
the results outlined in the introduction. First we calcu- 
late the two body correlations. 

III. TWO PARTICLE CORRELATIONS 

With expressions for the one and two particle density 
matrix we can easily calculate the two body correlation 
function, 

r^{ht2)~{h){t2) (22) 

where (ti) is the expectation value of the isospin U, of 
particle i, and so on. 

It may be helpful to emphasize that, since there is no 
momentum or spatial dependence in this model, the usual 
use of the correlation function to determine correlations 
between particles at different locations does not apply. F 
simply measures the correlation between isospins, but its 
time dependence is of interest. 

Expectation values may be obtained from density ma- 
trices in the usual way, e.g. 

{tit2) ^ Tt {pi2 \m1m2) {mim2\} , (23) 

or equivalently by first using the density matrix to obtain 
the probability P{mi,m2) that the isospin projections of 
the two particles are mi and 7712, and then calculating 

(^1^2)= "iim2P(mi,TO2). (24) 

mi ,m2 



Since the details of these calculations are given in the 
PhD thesis of Okuniewicz, we simply give the results 
here. 



A. The correlation function when both particles 
are initially in isospin-up states 

In the case where both particles are initially in the 
isospin-up state, Eq. (j20p applies. The final expression 
for the correlation function is rather long and not partic- 
ularly instructive. We give some results using it, allowing 
the interested reader to obtain the details from reference 

:2]. 

To analyse the correlation function we have plotted it 
as a function of the time for various numbers of isospin 
up particles, N, and isospin down particles, M. A subset 
of these simulations is shown in figure [T] and figure [2] for 
which iV = 101 and M is varied. 



r 




FIG. 1: A plot of the correlation function, F, using the initial 
conditon in which both particle 1 and particle 2 have isospin- 
up. In these plots iV = 101, M is varied such that N > M 
. The time is scaled as usual, r = gt{N + M). Note that 
the correlation function is positive for all cases. For the case 
= M the glitches in the plot are not real, rather they are 
an artifact of mathematica graphics 
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FIG. 2: A plot of the correlation function, F, using the initial 
condition in which both particle 1 and particle 2 have isospin- 
up. In these plots N = 101 and M is varied such that N < M 
. The time is scaled as usual, r — gt{N + M). Note that the 
correlation function is positive for all cases. 



The correlation function shows the same features as 
those found by Friedland, McKellar and Okuniewicz for 
the probability of one of the initial isospin up particles 
remaining in the isospin up state. Namely, 
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• When N ^ M, the system comes to equihbrium 
after some time. This features can be seen for 
the case N=101, M=101. Just as in the study of 
the equihbration of the single isospin by Friedland, 
McKellar and Okuniewicz, the equihbration time is 
incoherent. 

• When N ~ M is large the correlation function ex- 
hibits oscillations but the amplitude is small. 

• The correlation function is periodic on large time 
scales, with the same period as before. 

There are two new features to note: 

• the correlation is always positive, 

• when the equilibration time is clearly incoherent 
{N ^ M), the correlation is of order 10%, which is 
not particularly small. 

We defer discussion of these features until we have pre- 
sented the results for the case where the isospins are ini- 
tially anti-parallel. 
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FIG. 4: A plot of the correlation function, F, using the initial 
conditon in which particle 1 has isospin-up and particle 2 has 
isospin-down. In these plots A'^ = 101 and M is varied such 
that N < M . The time is scaled as usual, r = gt{N + M). 
Note that the correlation function is negative for all cases. 



Our study of the two body correlations in this solvable 
model show that, in precisely the case when the equilibra- 
tion occurs on an incoherent timescale, the correlations 
are not small. This leads us to ask the question: 



B. The correlation function when both particles 
are initially in anti-parallel isospin states 

In this case the density matrix of Eq. (PT|) is used. Fig- 
ure [3] and figure |4] show plots of the correlation function 
for N = 101 and various numbers of M. 
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FIG. 3: A plot of the correlation function, F, using the initial 
conditon in which particle 1 has isospin-up and particle 2 has 
isospin-down. In these plots N — 101, M is varied such that 
N > M . The time is scaled as usual, t = gt[N + M). Note 
that the correlation function is negative for all cases. For the 
case N — AI the glitches in the plot are not real, rather they 
are an artifact of mathematica graphics 

The features of this correlation function are similar 
to those found in the correlation function in the case of 
initially parallel isospins, and also by Friedland, McKellar 
and Okuniewicz for single particle probabilities. That 
is, the correlation function displays equilibration on an 
incoherent timescale, and also freeze-out and periodicity. 
In the case M ~ iV, the correlation is again of order 
10%. A difference from the parallel isospin case is that 
the correlation function is negative. 



Are the correlations big enough to induce 
significant corrections to the single particle 
kinetic equations, in spite of the incoherent 
time evolution? 

Since we are working with a solvable model, and thus we 
know the time evolution of the complete many particle 
density matrix, we can proceed to answer this question. 



IV. CORRELATIONS AND THE TIME 
DEVELOPMENT OF THE SINGLE PARTICLE 
DENSITY MATRIX 

A. The BBGKY Equation In the Model 

The time evolution of the density matrix, is described 
by the Von Neumann equation, 

= (25) 

We can split the Hamiltonian into its single particle and 
interaction constituents, 

^ 1 

^ = E^'+2E^^^' (26) 

i=l i=£j 

where K is the single particle Hamiltonian of the ith 
particle and Vij is the interaction Hamiltonian between 
the zth and jth particle. 

The rate of change of a single particle density matrix 
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is found by tracing over all the other particles, 



i 5] Tr, pg V 1 ^ Tr^y.i , ] 



[Xi,pW]+ETr,[yi„pi,] 

3 



(27) 

In the last line of Eq. (|27p we have employed the symme- 
try of our interaction Hamiltonian. 

Eq. ([27|) is the first quantum Bogolubov-Born-Green- 
Kirkwood-Yvon (BBGKY) equation ;30]. 

In the Friedman-McKellar-Okuniewicz model there is 
no single particle Hamiltonian. In this general discussion 
we retain such a Hamiltonian, which is responsible for the 
free space neutrino oscillations, but in numerical exam- 
ples it is omitted, although it can be introduced with just 
a little extra effort. We will frame this general discussion 
for the M flavor case. 

For the interaction Hamiltonian of our model it can 
readily be shown that the last term goes to zero, and 
the three particle density matrix does not contribute to 
the time dependence of the single particle density matrix. 
This is a significant simplification in the dynamics of the 
model. To demonstrate this represent the one, two and 
three particle operators on the flavor space on the basis of 
direct products of the extended Gell-Mann matrices A^, 
/i S {0, 1, . . . , A/"^ — 1}, which were introduced above. The 
superscript ['1 indicates that the Pauli martix or function 
refers to particle i, and similarly for multiple bracketed 
superscripts. Explicitly 





- -;cWa[^i 




(28) 








(29) 




2 " 




(30) 








(31) 








(32) 



Here the superscript refers to particle i, and we have 
employed a summation convention over repeated Greek 
indices. The few particle density matrices satisfy 



(33) 
(34) 



Note that Eq. ([34]) is valid for all iV — 1 possible values 
of j. From these follow the relationships 



' n — 



(35) 
(36) 
(37) 

pW = (38) 
The algebra of the extended Gell-Mann matrices is 

{A^,A,} = 2V„A« (39) 
[A^,A,] = 2i/^,,A« (40) 
XfjiXv — df^uKXK^iffji,vK,XK- (41) 

The coefficients d^^^, defined from the anticommutators 
are completely symmetric in their indicies, and the coeffi- 
cients f^vK, defined from the commutators are completely 
antisymmetric. These coefficients may, if required explic- 
itly, be calculated from 

1 —i 
d^PK = -^t: {{Xfi, Aj/} Xk) ffj.vK = ([A/ii Xy\ Xk) . 

(42) 

Using this algebra, the components of the final term in 
Eq. P?)) may be written as 



Tr,, 



\V - o'^^h - — vfelTrll^^l AW 



32 

lt,,(*W.AWA[^lAl^+AHAW(z/,,,^Agl) 



= 0, 



(43) 



showing that the BBGKY equation for p^^^ depends only 
on p(^^ and p^^^ . 



.dp[ 
' dt 



[ifi,pW]+5]Tr,[yy,pg)]. 



(44) 



In terms of the Gell-Mann matrix representation of the 
density matrices, Eq. is 



= /a. 



/3m 



(45) 



Note that, in the general case, Pg is constant, as it must 

be, but that 'Pl*! is not constant. Thus in general the one 
particle system is not a pure state, in that (p'^')^ 7^ p'^'. 
We have an explicit representation of the one and two 
body density matrices, in the case that Kj — 0, from the 
previous section, and it is a good check on our calcula- 
tions to verify that the Eq. (|44)) is indeed satisfied. This 
has been done, and as the calculation is rather involved, 
the interested reader is referred to the Honours thesis of 
JQ [3] for the details. 
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B. The Boltzmann equation and its failure in the 
model 



While it is reassuring that the exact BBGKY equation 
is satisfied in our exactly solvable model, we really need 
to test the analogue in the model of the neutrino kinetic 
equations. As shown in detail by Thomsonfo'l and McKel- 
lar and Thomson[ll| there are two key approximations 
in deriving the neutrino kinetic equations. 

• Two particle correlations can be neglected. 

• The time of a collision is negligible compared to the 
time between collisions. 

In the present model, as we have no spatial dependence 
in the interactions we cannot model the second approxi- 
mation. The first approximation can be implemented in 
the model by approximating the two body density matrix 
as the direct product of single particle density matrices, 



(2) 
Pi/ 



Pi ^ Pj 1 



(46) 



and substituting this approximation in Eq. (1441) . In this 
way we obtain a version of the quantum Boltzmann equa- 
tion 



/..J 4^1+ EOT 



(47) 



Before we specify the interaction appropriate to the 
model, we see that in the Boltzmann equation Eq. (IT71) 
the interactions now give a modified oscillatory be- 
haviour, as "precesses" about the vector with compo- 
nents /Cq -I- J2 



,^.^L1^AUndthus 



is a constant. 



which can be fixed to 1. With the general result that 
Vj/ = 1 this implies that {p^ 



l)^2 



P 



(1) 



In the Boltz- 



mann approximation the single particle density matrix 
represents a pure state. 

Now we note that in the model, the Gell-Mann matrix 
representation of the interaction is 

Vap = (t^aSaP with (/)o = 4g (TV + 1 ) /TV 

and (f)a = 2g for a 7^ 0, (48) 
so the model Boltzmann equation is 



where the restriction j i is no longer necessary as the 
contribution from j — i vanishes. From this it follows 
that T'I'l "precesses" freely, i.e. undergoes free os- 
cillations, and for the case studied above in which /C'*! 
vanishes, this vector is a constant. 



Further, if we specify the initial state so that all par- 
ticles at i = are in flavor eigenstates, as was done for 
TV = 2 above. 



■pW(t = 0) = unless a = n^-l, 
for 71 e {1, . . . ,TV}. 



Thus at t = 



dt' P 



0, 



and so at the shghtly later time St, 

(t = St)^0 unless a = n^ -1. 



Thus 



—Va =0 at all times. 
at 



(50) 



(51) 



(52) 



(53) 



This is precisely the result obtained by explicit 
calculation (3| in the case Af = 2, which is in complete 
contradiction with the results of Friedland, McKellar and 
Okuniewicz[l[, in which the exact non-trivial time devel- 
opment of p^^^ was computed. The Boltzmann equation 
is not satisfied in the model in the most dramatic way 
possible — all of the time dependence of p'^' is 
generated by the correlations — even though the 
correlation function, as we saw, has a magnitude of at 
most about 10%. 

Clearly the Boltzmann equation is not even approxi- 
mately valid in this model. 



V. CONCLUSION 

It is tempting to suggest that the breakdown of the 
Boltzmann equation in this exactly solvable model is 
closely connected to the structure of the model, and is 
not a general effect. In part the result is a consequence 
of the initial conditions. Even for Af — 2, ii we prepare 
some of the particles initially in an eigenstate of, say, t^ or 
the mass operator, and others (the flavor eigenstates) in 
eigenstates of tz , then the rhs of the Boltzmann equation 
is non- vanishing. Without further calculations it is not 
possible to say how well the Boltzmann equation is sat- 
isfied, but at least it is not completely and dramatically 
wrong. The way in which one can obtain a non- vanishing 
time dependence of the one particle density matrix from 
the Boltzmann equation is to have some of the initial 
neutrinos as one class of eigenstates, and the rest as an- 
other class of eigenstates. It is not easy to imagine how 
nature may have prepared such a state. 

The result is specific to the model interaction, and that 
interaction is not particularly realistic. An outstanding 
challenge is to generalize the interaction, maintaing the 
solvability of the model, and testing the Boltzmann equa- 
tion in this more appropriate model. 
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Regrettably, we have found that, while a coherent time Acknowledgements 
scale may indicate the importance of correlations, an in- 
coherent timescale is not a reliable indicator of the ab- 
sence of significant correlations. As a consequence this 
particular solvable model is not a good guide to the va- 
lidity or otherwise of the neutrino kinetic equations in This work has been supported in part by the Aus- 
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